Possible Z 2 phase and spin-charge separation 
in electron-doped cuprate superconductors 
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The SU(2) slave-boson mean-field theory for the tt'J model is analyzed. The role of next-nearest- 
neighbor hopping t' on the phase diagram is studied. We find a pseudogap phase in hole-doped 
(HD) materials (where t' < 0). The pseudogap phase is a (7(1) spin liquid (the staggered-flux 
phase) with a (7(1) gauge interaction and no fractionalization. This agrees with experiments on HD 
samples. The same calculation also indicates that a positive t' favors a Z2 state with true spin- 
charge separation. The Z2 state that exists when t' > 0.5J can be a candidate for the pseudogap 
phase of electron-doped (ED) cuprates (if such a phase exists). The experimental situation in ED 
materials is also addressed. 
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I. INTRODUCTION 

A complete understanding of the behavior displayed 
by high temperature superconducting cuprates is still 
lacking. In these materials, d-wave superconducting 
(dSC) samples are obtained upon doping the parent com- 
pounds. The latter are Heisenberg antiferromagnets. 
The presence of low energy properties consistent with 
standard BCS theory in dSC samples is also generally 
undisputed. However, the phase-space region interven- 
ing between the two aforementioned regimes hosts an un- 
usual phenomenology - and a quite debated one as well. 
In this paper, we will adopt the point of view that the 
unconventional behavior in underdoped cuprates is evi- 
dence of doped Mott insulator physics. To capture the 
important properties of these materials, which are then 
controlled by the large Mott gap and the two-dimensional 
spin interactions in the copper-oxide planesi, we will use 
the slave-boson approach developed in Refs. and0. 

In underdoped electron-doped (ED) materials the an- 
tiferromagnetic (AF) phase is very robust. However that 
is not the case in their hole-doped (HD) counterparts. 
Instead, a pseudogap metallic normal state appears in 
underdoped samples. This is a paramagnetic state with a 
d x 2_ y 2 gaped spectral function. Reconciling the gap with 
the strong local AF fluctuations is not trivial though. 4 
One way to understand it is offered by the slave-boson ap- 
proach to the t J model. In particular, in Ref. |3|a transla- 
tion symmetric pseudogap metallic state (the staggered- 
flux or s-flux state) was proposed. It contains spinons 
and holons interacting via a (7(1) gauge field. This state 
was shown to bear enhanced AF fluctuations compati- 
ble with experiments* 4 " The (7(1) gauge fluctuations play 
a key role in the AF instability at very low doping as 
well as in the destruction of quasiparticle peaksA& How- 
ever, as discussed in Ref. 0, the appearance of the spin- 
pseudogap below a certain energy and the existence of 
well defined clcctronlikc quasiparticles in the nodal di- 
rection above T c suggest that spin and charge recombinc 
before the s-flux state reaches the dSC state. 

There is a competing scenario to describe this crossover 
from the s-flux to the dSC phase which involves the emer- 



gence of an intervening state with a Z2 gauge structured 
This state would result from the condensation of spinon 
bilinearsiSi2ii£ The dSC state would then follow from 
holon condensation. Z2 spin liquids are very interesting 
states with true spin-charge separation. They were pre- 
dicted theoretically over ten years agoAS Yet they lack 
experimental realization. The present work proposes to 
study the conditions favoring the emergence of these spin 
liquid states as well as to discuss their relevance for the 
cuprate superconductors. 



Despite the effort to check the applicability of Z2 spin 
liquids in the context of high T c superconductivity, so 
far all the experimental tests gave negative resultsiii* 1 ^ 
However, this does not imply that Z2 states do not exist 
in high T c samples. In this paper we find that Z2 states 
are indeed unlikely to exist in HD compounds. But a Z2- 
linear state {i.e. with linear gapless spinon excitations) 
is likely to appear as the pseudogap metallic state for 
ED materials. In superconducting (SC) samples without 
AF order spin liquid pseudogap signatures may emerge 
once SC coherence is destroyed by thermal fluctuations 
or an external magnetic field. In the first case such a 
behavior is still to be reported. Experimental evidence 
in the magnetic field driven normal state, however, seems 
to be consistent with Z 2 pseudogap phenomenology! 13 i 14 



The t J model has a particle- hole symmetry. In order to 
break it we introduce the next-nearest-neighbor (NNN) 
hopping t' term. Taking both signs of t' covers both ED 
(f > 0) and HD (t' < 0) cases. In Sec. |H]we introduce 
the SU(2) MF theory for the tt'J model. In Sec. UlTlwe 
present the resulting MF phase diagrams. For t' < our 
calculations are consistent with the previous studies for 
t' = where only the s-flux state was obtained. How- 
ever we find that NNN hopping can stabilize a ^-linear 
state in the tt'J model for values of t' > 0.5 J. These 
results show that unfrustrated hopping favors fractional- 
ized phases (Sec. II V|) . In Sec. |V|we establish the link 
between our results and some properties of ED cuprates. 
We draw our main conclusions in Sec. IVII 



2 



II. 517(2) MF THEORY FOR THE tt'J MODEL 



The tt'J model is given by the Hamiltonian 
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where sums are taken over pairs of sites (ij). The ex- 
change coupling Jij is J for (ij) nearest-neighbor (NN) 
sites while the hopping parameter tij equals t for (ij) 
NN sites, i' for (ij) NNN sites. 5; is the spin on site 
i, c\ and rii = cjcj are electron creation and occupation 
number operators and P projects out doubly occupied 
sites. 

To deal with the no double occupancy constraint 
the electron operator is decoupled according to a-\ = 



Si = 



44) 



and 



+ *fe4) and C U ~ 73V 
the spin operators are expressed as Si = \j\ a T a j3fi(} 
where the SU{2) doublets *j = [ipj^] = [flxfa] and 
™i = [44] represent spinons (fermionic particles carry- 
ing spin 2) and holons (bosonic particles carrying unit 
charge) respectively and f are the Pauli matrices. These 
operators were introduced in Ref. to keep the SU(2) 
gauge structure of H U 'j even away from half-filling. Par- 
ton operators introduce unphysical degrees of freedom 
which are projected out by implementing the constraint 
(i.e. physical states are SU(2) 



hjfhi 



singlets on every site) . Introducing a Lagrange multiplier 
per lattice site (aoi) the problem is reduced to one of lat- 
tice gauge theory^ Performing a Hartree-Fock type of 
decoupling of Hh f and Hhop (neglecting H' for the time 
being) and implementing the projection constraint at MF 
level only, the following MF Hamiltonian is obtained: 
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where fib is the chemical potential of holons. 

Despite ignoring gauge fluctuations, we expect MF the- 
ory to account for some qualitatively right features. In- 
deed, we have the following: 



i) The energetics of the model is correctly captured at 
MF level. This fact is suggested by the SU(2) slave- 
boson MF theory phase diagram for the t J model - it in- 
cludes the pseudogap regime (s-flux state) that turns into 
the dSC state as holons become coherent, as well as the 
strange-metal and Fermi-liquid regimes^ In particular, 
the pseudogap metallic state was predicted by the slave- 
boson approach prior to experimental observation^ 

ii) The nature of fluctuations in each phase is also con- 
sistent with experiments. For instance, the U(l) gauge 
structure in the s-flux state accounts for the spin excita- 
tion spectrum in underdoped HD materials^ It also leads 
to the lack of well defined quasiparticles. 5 In the strange- 
metal regime SU (2) gauge fluctuations are responsible 
for the observed non-Fermi liquid behavior; 

Hi) MF states have properties that survive projection 
even in the presence of long range U(l) gauge fields. 
Take the s-flux state, which includes both staggered cur- 
rent and d-wave pairing fluctuations, as an example. Re- 
markably, in Ref. ^] staggered vorticity correlations were 
reported to emerge from Gutzwiller projecting the dSC 
wave function. There is also numerical evidence that Zi 
MF states lead to fractionalized phases after performing 
Gutzwiller projection. 17 



III. 



MF PHASE DIAGRAMS 



To implement the self-consistent MF equations for the 
tt'J model, the uniform ansatz considered in Ref. 0j 
must be extended to include the NNN MF parameters 
Ui,i+x±y- Here we only consider translation invariant 
ansatse which do not break any symmetry. In our MF 
calculation we use the following ansatz which properly 
comprises both the s-flux and dSC states: 



Ui,i+$ = xt 

U^+y = XT 3 

Ui,i+£±y = T 7 " 3 ! 

a 3 , ± 0. 



TjT 1 



(4) 



Depending on the values of \> Vi an d 7, the above ansatz 
can describe phases with the same symmetry but differ- 
ent quantum orders i& These phases may contain SU(2), 
U(l), or Zi gauge fluctuations. Some phases display a 
large Fermi surface while others have only Fermi points. 
These different MF states and their labels are summa- 
rized in Table HJ 

The MF phase diagrams are determined for t = 3 J 
and t' between — 1.5J and 1.5J. These values are rep- 
resentative of the cuprates (note that t' is close to J for 
ED samples). Holons may condense at nonzero temper- 
atures due to interactions or to the projection constraint 
(as proposed in Ref. 0). However, at MF level holons are 
free and conventional holon condensation only occurs at 
T = 0. In our calculations we take T > and, hence, 
only consider the unconventional mechanism. 
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TABLE I: Labels for MF phases, p is the density of condensed 
holons. All other parameters are defined in the main text 
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FIG. 1: (a) (t',x) phase diagram for H M f at T = + . (b) 
(t',x) phase diagram at T — + after including H' and JJ" 



The resulting (t 1 , x) phase diagram for T = + is shown 
in Fig^a). For t' < 0.5 J the s-flux state survives all the 
way to T — + - it becomes the dSC state for T = 0. 
However, for t' > 0.5J new MF states are obtained. In 
the underdoped regime, where the SU (2) theory was pro- 
posed to be relevant, the Z2 phase emerges. When T = 
holons condense in the band bottom and the Z2 state 
changes into the dSC phase. For high enough doping we 
find the dSC state even at nonzero temperature - result- 
ing from projection constraint driven condensation. As 
the pseudogap closes (77 — > 0) a state with U(l) gauge 
structure and gapless spinons arises (Ul). Upon holon 
condensation it becomes a Fermi liquid (FL). 

For t' < 0.5 J we obtain the (x,T) phase diagram re- 
ported previously^ The MF finite temperature phase di- 
agrams are illustrated for t' = 0.8 J and t' — 1.2 J in 
Figs. Efa) and [2Jb) respectively. As, by construction, 
the AF phase is absent, the ^-linear state appears as 
the dominant phase at low doping. ED samples can have 
three different normal metallic states: (a) a Z2 state with 
a d-wave pseudogap, true spin-charge separation and a 
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FIG. 2: (x,T) MF phase diagrams for (a) t' = 0.8J and (b) 
t' — 1.2J. In (c) and (d) qualitative phase diagrams, based on 
the MF results from (a) and (b) respectively, are proposed. 
These include the experimentally observed AF phase and fi- 
nite temperature holon condensed phase due to holon inter- 
action. The dashed line describes the transition into the Z2 
state in case it is not completely absorbed by the AF phase. 
Above T c the Z2 state describes pseudogap behavior. In (d) 
this transition may describe the reported transition from dSC 
to an extended s-wave SC (sSC) phase in P^-zCe^CuO^-,, 
and La2-xCe :!; Cu04-j / (Ref. Il8f) . 



nontrivial topological order(2*S^ (b) a strange metal (the 
Ul phase) with a large Fermi surface and a U(l) gauge 
interaction, and (c) a FL phase. 

The above discussion of the MF phase diagrams is 
qualitatively correct even after including the dropped H' 
term. In terms of spinons and holons, H' has the form 

J 



2 

J 
16 



E E hthihl^h+z + 0(x 3 ), (5) 



t u—x,y 



(which is oc x 2 at the MF level). While the second term 
correctly describes the attraction between holons result- 
ing from the background with AF fluctuations, the first 
one should have the opposite sign to correctly account 
for the attraction between spinons and holons caused 
by gauge fluctuations. Hence, H' should be considered 
together with an effective contribution arising from the 
gauge interaction. At MF level the constraint is imple- 
mented only globally, i.e. \^2iQi) = 0- Including the 
following effective interaction, H" = £\ . (Qj — 

Qi + a) 2 , the fluctuations of (Qi) 2 away from zero are 
reduced. The prefactor in H" is such that the ex- 
change constant is renormalized to J e ff — 1-5 J (note 
that Gutzwiller projection is known to increase the effec- 
tive exchange constant) . H" also makes holons more 
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massive. Including H' and H" at the MF level extends 
the region where the Z2 state emerges (Fig^b)). Hence 
gauge fluctuations stabilize the Z2 phase. As expected, 
projection has a smaller effect in this phase. 

The above results support that the Z 2 -linear spin liq- 
uid introduced in Sec. [n]is a competitive state for higher 
values of t' and low doping levels. In particular, the 
lower bound t' > 0.5 J makes this state relevant for ED 
cuprates. 



IV. SPIN-CHARGE SEPARATION: THE ROLE 
OF 

The SU(2) slave-boson approach is designed to address 
the low doping regime. As it is clear from Fig^ for low 
doping there are two dominant phases, namely, the s-flux 
state and the Z2 state. The former becomes the latter 
as t' is increased above a certain critical value. This pa- 
rameter driven transition links two states with the same 
symmetries (they both describe totally symmetric spin 
liquids). However they have different quantum orders^S 
This statement can be made more explicit by emphasiz- 
ing the quite different nature of their low-energy effective 
behavior. 

In the s-flux state massless Dirac fermions (spinons) 
and charged bosons (holons) interact with a long-range 
U(l) gauge field. The gapless nature of these interacting 
excitations is protected by the quantum order>iS Still, 
the gauge field has a drastic effect on the nature of the 
excitations - this state is in fact a manifestation of a 
two-dimensional Luttinger liquidii^ which has no well 
defined quasiparticles. 

On the other hand, in the Z2 state the gauge interac- 
tion becomes short-ranged. Hence, the (linearly dispers- 
ing) spinons and holons are well defined quasiparticles. 
For that reason there is true spin-charge separation. The 
remnant discrete Zi gauge structure reflects the topolog- 
ical order of the stated This topological order is related 
to the ground-state degeneracy of the system when it is 
embedded in a manifold with nontrivial topology. Such 
degenerate ground-states are locally similar but have dif- 
ferent global (topological) properties. 

The different gauge structure in the two states (U(l) 
and Z2) naturally gives rise to qualitatively distinct spec- 
tra for the collective modes. Indeed, the transition is ac- 
companied by the opening of a gap in the fluctuations 
around the MF saddle-point^ 

Physically, the transition from the s-flux state to the 
Z2 state corresponds to the emergence of coherent di- 
agonal (intrasublattice) charge carrier hopping - i.e. 
(cfej+fi + c|, fi Cj), with u — ±x ± y, becomes nonzero. It 
results from a combination of two factors: (a) increasing 
t' unfrustrates diagonal hopping and (b) electron/hole 
intrasublattice hopping is not frustrated and does not 
frustrate the background AF correlations (unlike inter- 
sublattice hopping). These play a crucial role in the dy- 
namics of charge carriers (and of spin degrees of freedom 



as well). 

In the s-flux state intrasublattice hopping is depleted 
by the negative t'. As charge carriers hop between dif- 
ferent sublattices they interact strongly with the back- 
ground AF fluctuations. Hence charge carrier motion 
and AF correlations frustrate each other. That is the 
motive underlying the staggered current and staggered 
chiral spin correlations The presence of such correla- 
tions in the context of the t J model was reported in exact 
diagonalization studies^ Moreover, charge carrier den- 
sity and vorticity correlations are proportional to each 
other. 16 - 22 The emerging picture consists of charge carri- 
ers of opposite vorticity attracting each other in a back- 
ground of staggered chiral spin fluctuations. This attrac- 
tion eventually induces dSC coherence. 

In the Z2 state charge carriers hop coherently along 
the diagonals and the staggered current correlations de- 
crease. In particular, the energy difference between the 
dSC and staggered current states is seen to increase with 
t'. Therefore, staggered current fluctuations are reduced 
in the Z2 state. Exact diagonalization studies in doped 
tt' J ladders show the same trend - i.e. staggered currents 
are stabilized by decreasing t'm& According to our calcu- 
lations, spinon pairing (77) is also reduced by the emer- 
gence of coherent diagonal hopping (therefore, charge 
carrier pairing may be expected to decrease as well) . This 
is further consistent with the fact that the ^-vortex gap 
is physically connected to the reduction of staggered chi- 
ral spin fluctuationsi2fl*2i 

Topological order in the Z2 state is associated with 
the intrasublattice hopping parameter 7. Studies of 
Gutzwiller-projected wave functions in Ref. fur- 
ther corroborate the importance of the NNN hopping 
term (i.e. a non- vanishing 7) in realizing true fraction- 
alization. Consequently, for underdoped samples, the 
^2-vortex gap in this fractionalized metallic pseudogap 
phase is set by 7 which is oc x. This is qualitatively dif- 
ferent from the suggestion of Senthil and Fisher where 
the ^-vortex gap is set by the pseudogap scale (rj)H 

The arguments above support and delineate how the 
parameter t' induces spin-charge separation. However, to 
fit experimental band structure, an extra hopping term 
(i") must be included^ This requires an extension of 
the ansatse considered in our analysis. Nevertheless, in a 
tight-binding model with t > 0, t' > and t" < (as it is 
the case for ED cuprates), t" is seen to frustrate hopping 
along the lattice bonds while unfrustrating hopping along 
the diagonals - in this sense enhancing the effect of t'. 



V. APPLICATION TO ED CUPRATES 

As mentioned before, in ED cuprates the AF state 
is very stable. Such stability can be addressed within 
the tt' J model framework. Indeed, an increasing t' un- 
frustrates intrasublattice hopping and renders the charge 
carriers less effective in frustrating the background AF 
correlations^ 28 However, our MF results suggest that, 
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if AF order can be destroyed by tuning parameters in 
the Hamiltonian, the resulting spin liquid may contain a 
pseudogap and correspond to the Z2 state. 

A distinctive experimental feature of the Z2 phase is 
the presence of a d-wave pseudogapi^L Well, in the differ- 
ent Nd 2 -xCe x Cu0 4 -y, PT 2 - x Ce x Cu0 4 - y (PCCO), and 
La 2 _ 3; Ce :r Cu04- a (LCCO) SC thin films, the SC gap 
(~ 5— lOmeV) is reported to survive in the magnetic field 
driven normal stated The value of the normal state gap 
is comparable to the SC gap and decreases with doping. 
Moreover, the low temperature normal state in PCCO SC 
crystals was reported to violate the Wiedemann-Franz 
law^ This provides evidence for a non-Fermi liquid state 
compatible with spin-charge separation. According to 
our MF results, the aforementioned normal state may 
realize the Z2 state. On the other hand, if the Z2 state 
is to exist above T c it probably only covers a small frac- 
tion of the phase diagram due to the robust AF phase 
(Figs Etc) andEfd))^ 

Despite still disputed, there are a number of experi- 
ments reporting dSC order in ED samples^ In our pic- 
ture, the d x 2_ y 2 gap in the SC phase is inherited from 
spinon pairing in the Z2 state. Recently, however, a SC 
pairing symmetry transition across optimal doping, from 
d-wave to extended s-wave, was reported in both LCCO 
and PCC04& We suggest the extended s-wave pairing 
may result in the FL state (Fig|2Jd)) from exchange of 
collective modes - as the interaction between the Fermi 
sea and the massive gauge-fields raises an instability in 
the s-wave pairing channel^ We remark that this mech- 



anism is unrelated to the one underlying dSC. 



VI. CONCLUSIONS 

In this paper we use a self-consistent SU (2) slave-boson 
approach to study the role of NNN hopping in the context 
of the tt'J model. We find that t' induces a transition 
between the s-flux state (an algebraic spin liquid^) and 
the Z2 state (a spin liquid with true spin-charge sepa- 
ration). This transition is interpreted as the emergence 
of coherent intrasublattice hopping in the Z2 state. The 
parameter t' is also found to decrease staggered current 
and staggered chiral spin fluctuations, as well as to re- 
duce pairing between charge carriers. 

The range of parameters involved leads us to propose 
that a pseudogap metallic state in ED cuprates is likely 
to be described by a Z2-linear phase. In contrast, the 
pseudogap state for HD samples is found to be described 
only by the s-flux phase with U(l) gauge interaction. 

Realizing fractionalized states in real systems is a ma- 
jor challenge. The present work suggests that the search 
for a pseudogap metallic state in ED cuprate samples 
may achieve such a goal. 
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